Estimating Total Solar Irradiance during the 21st century 



Victor Manuel Velasco Herrera,^ Blanca Mendoza^ and Graciela Velasco Herrera.^ 

^ Departamento de Ciencias Espaciales, Instituto de Geofisica, Universidad Nacional Autdnoma de Mexico, 
Ciudad Universitaria, Coyoacdn, 04510, Mexico D.F., MEXICO and 
^Centro de Ciencias Aplicadas y Desarrollo Tecnologico, Universidad Nacional Autdnoma de Mexico 

The reconstruction and prediction of solar activity is one of the current problems in dynamo theory 
and global climate modeling. We estimate the Total Solar Irradiance for the next hundred years based 
on the Least Square Support Vector Machine. We found that the next secular solar minimum will occur 
between the years 2003 and 2063 with an average of 1365.4 W/w? close to the Dalton or Modern minima. 
We calculate the radiative forcing between the modern maximum and the 21^^ century minimum to be 
-id.lW/m^. 



The modern comprehension of solar variability possibly 
began when Schwabe published the periodicity of sunspots 
in 1843 [IJ. In 1894, Maunder published a discovery 
that has maintained the Solar Physics in an impasse. In 
his famous work on "A Prolonged Sunspot Minimum" 
Maunder wrote [2J: 

"The sequence of maximum and minimum has, in 
fact, been unfailing during the present century [..] and yet 
there [..], the ordinary solar cycle was once interrupted, 
and one long period of almost unbroken quiescence 
prevailed" 

The Maunder's discovery went unnoticed and forgotten 
until 1976 when Eddy brought it to light again [3|. The 
existence of prolonged solar minima has been one of the 
most controversial questions in Solar Physics. However 
the possibility of prediction of new periods of diminished 
solar activity is even more controversial. 

To estimate future solar activity, several methods have 
been used, for instance dynamo models, spectral meth- 
ods, regression methods or neuronal network methods H. 
These estimations have in common that they are applied to 
short-time reconstructed series and that they do not discuss 
the relative accuracy of the methods. 

The fluctuations of the solar time series is a tool that 
helps to study the solar magnetic field as well as to un- 
derstand the solar dynamo. These fluctuations can occur 
for instance in the amplitude, phase, frequency, energy and 
power of the solar phenomena. 

The majority of solar activity analysis focuses on fluc- 
tuations of the amplitude. In this paper, we propose to 
consider not only the fluctuations in amplitude but also in 
the power of the Total Solar Irradiance (TSI) as a physi- 
cal measure of the energy released by the solar dynamo, 
which contributes to understand the nature of 'profound so- 
lar magnetic field in calm". Regardless of the mechanism 
that produces solar activity minima (stochastic, chaotic, in- 
termittent or quasi-periodical processes), the study of these 
minima is very important for the solar dynamo theory, 
as well as for its impact on solar-terrestrial relationships 



Recent studies suggest that the mid- term (1 — 2 years) 
and the secular periodicities are the product of chaotic 
quasi-periodic processes and not of stochastic processes or 
intermittent process ||71[8l. Different spectral analysis of 
solar activity series E [9l show several significant long- 
term periodicities. It is also known that the solar cycle 
(Schwabe periodicity [1 J) varies cyclically with a mean pe- 
riod of about 11 -years and the magnetic cycle (Hale cycle 
[ 10|) with a mean period of about of 22-years. This behav- 
ior motivates attempts to predict solar activity, especially 
now that an unexpectedly low activity solar cycle 23 oc- 
curred and could be the sign of the beginning of a new sec- 
ular solar minimum [8|. The behavior of the solar cycle 
23 minimum has shown an activity decline not previously 
seen in the past cycles for which spatial observations exist 

EHIl. 

The descending phase and minimum measurements of 
solar cycle 23 show that in particular the TSI has fallen 
below the previous two solar minima values: the mean 
PMOD composite TSI for September 2008 is 1365.26 ± 
0.16 Wm-^ , compared to 1365.45 Wm-^ in 1996 or 
1365.57 Wm-^ in 1986 113. 

Over the 11-years solar cycle, TSI variations of ^0.07% 
have been observed between solar minimum and maxi- 
mum |[T6ll . This modulation is mainly due to the inter- 
play between dark sunspots and bright faculae and network 
elements! 17]. Studies using cosmogenic isotope data and 
sunspot data ^ [TSl indicate that currently we are within a 
grand activity maximum which began after ^1930. 

Studying the solar wind, the interplanetary magnetic 
field strength and the open solar flux over the past cen- 
tury, Lockwood [19 ] found that all three parameters show 
a long-term rise peaking around 1955 and 1986 and then 
decline, yielding predictions that the grand maximum will 
end in the years 2013, 2014 or 2027 depending on the pa- 
rameter used. 

Other works indicate that the current maximum will not 
last longer than two or three solar cycles more ll20l . Fur- 
thermore, it has been suggested that a Dalton-type mini- 
mum has already began in solar cycle 23 reaching 1 211 1221 
to solar cycles 24 and 25, while the Solar Cycle 24 Pre- 
diction Panel indicates a lower limit of 90 =b 10 for the 
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maximum sunspot number f23] of solar cycle 24. 

To project the TSI for the next hundred years, we use 
a method based on Least Squares Support Vector Ma- 
chines (LS-SVM) with Nonlinear Autoregressive Exoge- 
nous (NARX) model and with radial basis function (RBF) 
kernel that allows a better precision in the estimation of the 
future values of a time series |[24l[25ll . 

Since the instrumental TSI records began in 1979, we 
work with reconstructed TSI series. Long-term reconstruc- 
tions of TSI EI show epochs of maxima and minima 
when substantial changes in the TSI occur. These changes 
can contribute to climate variability [ 28 1 and new estimate 
of the TSI for the 21st century, will have an important im- 
pact on climate modeling. 

We tested three TSI recent reconstructions. Two of them 
use a model based on the evolution of the Suns total and 
open magnetic flux: The Wang reconstruction [29| con- 
sider differential rotation, supergranular convection and 
meridional flow; this model is used to derive two TSI re- 
constructions, one with and one without a secularly vary- 
ing background; from 1850 onwards these reconstructions 
are the recommended solar forcings for the fifth Coupled 
Model Intercomparison Project 20th century simulations 
[301. 

The Krivova reconstruction |[3H relies on time con- 
stants representing the decay and conversion of the dif- 
ferent photospheric magnetic flux components. The Sten- 
hilber reconstruction! 26], uses the observationally derived 
relationship between TSI and the open solar magnetic field, 
the latter obtained from the cosmogenic radionuclide ^^Be. 

However, the best reconstruction to apply our method is 
the Krivova reconstruction [ 31 1 because it has a better tem- 
poral and spectral resolution compared to the Stenhilber 
reconstruction [26 ] and the correlation between the recon- 
structed and the PMOD composite is better than the Wang 
reconstruction ll29l . 

From 1610 to 1975 we use the Krivova reconstruction 
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FIG. 1. Time series, (a) The TSI-KRIVOVA-PMOD (red line) 
from 1610 to 2010 superposed on the times series estimated with 
the Least Squares Support Vector Machines (LS-SVM) between 
1610 and 2010 (black line), (b) The LS-SVM model time series 
(black line) between 1610 and 2100 
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FIG. 2. (a) Upper panel, TSI-KRIVOVA-PMOD; right panel 
spectral wavelet analysis; left panel global wavelet, (b) The same 
as in (a) but for the LS-SVM model. 



ySTj and from 1976 to 2010 we use the PMOD composite. 
In this work we use the TSI-KRIVOVA-PMOD between 
1610 and 2010. 

To search for the future TSI values, the NARX LS- 
SVM was trained with 80% of random data of the TSI- 
KRIVOVA-PMOD, obtaining a mean squared error (MSE) 
of 0.0084, the testing of the remaining 20% presents an 
MSE of 0.0073. 

We plot the TSI-KRIVOVA-PMOD with red line and 
the NARX LS-SVM model with black line in Fig. la. 
It is clear that the NARX LS-SVM model (black line) 
reproduces very well the TSI-KRIVOVA-PMOD; the lin- 
ear correlation coefficient between the TSI-KRIVOVA and 
the NARX LS-SVM model series for 1610-1975 is r = 
0.9969, while r value for the PMOD and the NARX LS- 
SVM model for 1976-2010 interval is 0.9959. 

In Fig. 16 we present the NARX LS-SVM future TSI 
estimation (black line) between 2011 and 2100 with stan- 
dard deviation a = 0.4318 Wm~'^. We notice a de- 
creasing trend of the TSI between 2003 and 2063, co- 
inciding with other types of prediction adopting different 
methods|[8l[T9l[20l|22J. 

However, it is not enough to conclude that the amplitudes 
among the reconstructed, the composite and the modeled 
time series are similar, it is also necessary to compare the 
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spectral characteristics. 

We apply the wavelet analysis using the Morlet function 
[321 to quantify the TSI time series till 2100 A.D. and 
to analyse local variations of multiple periodicities. This 
method provides a higher resolution of periodicity, allows 
us to calculate the phase, and to filter the TSI in bandwidths 
L331. 

Also, to calculate the confidence level we used the nor- 
malizepdf function, in this way the TSI will have a gaussian 
distribution! 34 1 . Wavelet meaningful periodicities (confi- 
dence level greater than 95%) must be inside the cone of 
influence (COI), which is the region of the wavelet spec- 
trum outside which the edge effects become important [ 35 1. 
We also include the global spectra in the wavelet plots to 
show the power contribution of each periodicity inside the 
COI[36|. 

We established our significance levels in the global 
wavelet spectra with a simple red noise model (increasing 
power with decreasing frequency [37 1). We only took into 
account those periodicities above the red noise level. 

In Fig. 2 we show the wavelet analysis for the TSI- 
KRIVOVA-PMOD (Fig. 2a) and the NARX LS-SVM 
model (Fig. 26). In the central panel of Fig. 2a the 11- 
years cycle appears attenuated during the secular minima 
and is stronger since ^1875, while the ^120-years cycle 
keeps a more or less uniform power along all the time-span 
considered, both periodicities appear above the red noise 
level in the global spectra (right panel). Fig. 2b shows 
roughly the same spectral evolution of Fig. 2a but the two 
peaks above the red noise level are the 11 -years and the 
240-years. Also Fig. 2b indicates that the 11-years cy- 
cle will be attenuated during the next 60 years (to 2063), 
which is a characteristic of a secular minima. 

In Fig. 3 we show a further analysis of the amplitudes 
(black areas) and phases (blue lines), that allow us to quan- 
tify the starting and ending of a cycle of the NARX LS- 
SVM model and were obtained using the inverse wavelet 
transform [32J. Fig. 36 shows the 11-years cycle, as 
sunspots were so scarce during the "Prolonged Sunspot 
Minimum", it has been suggested that the solar dynamo 
stopped [5, 6|. The figure clearly shows the 11-years am- 
plitude attenuated during the Maunder, Dalton and Modern 
minima (negative phase of the periodicity of 120-years in 
Fig. 3c). The amplitude of the solar cycle never disappears 
completely, this is particularly so during the Maunder min- 
imum [38|. After the year 2000 the peak amplitude of the 
cycle tends to decrease, in fact during the 21st minimum 
(2003—2063) it is similar to the Dalton or Modern minima 
(1883-1940). 

The phase also shows an amplitude modulation, it 
presents an inversion between ^1780 and 1790, as this 
in not a particularly maximum or minimum time, we sug- 
gest that it comes probably from a shortcoming of the TS- 
KRIVOVA reconstruction in this time-interval. 

From 2011 to 2100, the phase does not show any other 
inversion, probably indicating the good quality of our es- 
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FIG. 3. Phases (blue lines) and amplitudes (black areas) associ- 
ated to the main periodicities of the LS-SVM model. The light- 
gray shading indicates the time-span of the secular minima, (a) 
Modeled time series . (b) The 11-yrs cycle, (c) The 120-yrs cycle, 
(d) The 240-yrs cycle. 



timation. Fig. 3c shows the 120-years cycle, we notice 
that it is the negative phase of this periodicity that coin- 
cides with the secular minima: Maunder between the years 
1637-1704, Dalton between the years 1771-1829 , Mod- 
ern between the years 1883—1940 and the 21st century 
minima between the years 2003 and 2063. 

Based on this periodicity we calculate the average TSI 
for the secular minima: Maunder 1364.8 Wm~^, Dalton 
1365.2 Wm-^ Modern 1365.5 Wm-^ and the 21 st cen- 
tury 1365.4 Wm-^ ± la = 0.4318 Wtr-'^, this TSI 
value is between the Dalton and Modern minima. Again 
we observe that the phase does not change after 2010. 

The Modern maximum noticed between the years 1940 
and 2003 has an average TSI of 1366 Wtr-^, and the 21st 
century minimum has a TSI average of 1365.4 Wm~'^, 
this implies a negative radiative forcing of ^—0.1 Wm~'^. 
The radiative forcing for the Maunder minimum to the 
Modern maximum[31 1 is ^0.21 Wm~'^, then the radiative 
forcing associated to the 21st century minimum is almost 
half of that forcing. 

Finally, Fig. 3d shows the 240-years cycle, it seem that 
according to the lag between the 120-years and 240-years 
cycles we have different amplitudes of the secular solar 
minima, for instance, when its negative phase coincides 
with the 120-years cycle we have the deepest minima, like 
the Maunder minimum. Again, the phase of this periodic- 
ity holds for the next 100 years. 

To decide when the solar activity is "high" or "low", we 
calculate the power of the TSI as a direct indicator of en- 
ergy released by the solar dynamo and the level of activity 
for each solar cycle. We use the mean power value of the 
PMOD composite (1976 - 2010) to calculate the anoma- 
lies for each cycle. The power anomalies are normalized 
and appear as blue bars in Fig. 4, which also presents the 
TSI (solid line). There are positive power anomalies during 
solar cycles 21 and 22, coinciding with the end of the Mod- 
ern maximum. Between solar cycles 23 and 30 the power 
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FIG. 4. We use the mean power PMOD composite (1976-2010) 
to calculate the power anomaly for each cycle. The power anoma- 
lies are normalized and appear as blue bars. The TSI is the 
solid line. Negative (or Positive) power anomalies is 'Low" (or 
"High") solar activity 



anomalies are negative, coinciding with negative phase of 
the 120-years periodicity. This reinforces the result of Fig. 
3c suggesting that this periodicity is closely associated with 
secular minima. According to the power of the anomalies, 
solar cycles 23-25 and 30 could present lower activity than 
cycles 26 to 29, regardless of the fact that the peaks of cy- 
cles 27 and 28 are the lowest. 

The calculated power anomalies show that low solar 
secular activity occurs when there are negative anomalies 
and high solar secular activity appears with positive 
anomalies. It is possible that the zero in the anomalies, 
represents the normal state of the dynamo. The "Prolonged 
Sunspot Minimum" discovered by Maunder, represents a 
phase of solar history and corresponds to a special state 
of the dynamo when it is working well below its average 
power. 
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